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Abstract. We consider the vacuum energy for a configuration of two cylinders and 
obtain its asymptotic expansion if the radius of one of these cylinders becomes large 
while the radius of the other one and their separation are kept fixed. We calculate 
explicitly the next-to-leading order correction to the vacuum energy for the radius of 
the other cylinder becoming large or small. 
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1. Introduction 



q ■ During the past years remarkable progress was made in the field of Casimir effect pQ and 
in the calculation of the vacuum interaction energy and the Casimir force for separated 
bodies of nontrivial shape. It turned out that it is possible to write down a representation 



of this vacuum energy in terms of a functional determinant which does not contain 
ultraviolet divergences. The first applications had been made for the configuration of 
two parallel cylinders and two spheres [21 [3] . From here, the corresponding expressions 
for a cylinder or a sphere in front of a plane follow for symmetry reasons. An alternative 
approach by [I] should be mentioned which finally delivers the same results. 

The approach using functional determinants allows for a direct numerical evaluation 
since all involved summations and integrations do converge. Especially for large 
separation between the interacting bodies the convergence is rapid. It is also easy to 
obtain an asymptotic expansion in this limit. It corresponds to a dipole approximation 
where only the lowest orbital momenta are involved. In the opposite limit of small 
separation the situation is more complicated. Here arbitrarily high orbital momenta give 
a significant contribution and reliable numerical results are hard to obtain. However, 
the asymptotic expansion for small separation can be calculated analytically. For a 
cylinder in front of a plane this was done in [5] and for a sphere in [6]. As a result, 
for these configurations, the Proximity Force Approximation (PFA) was re-confirmed 
and the first correction beyond was calculated. These results were partly confirmed by 
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Figure 1. The configuration of two cylinders 



numerical approaches. First we mention the remarkable world line method [3, [8] , which 
confirmed the above mentioned results for the case of Dirichlet boundary conditions. 
These results were also confirmed by extrapolation of the numerical evaluation of the 
functional determinant from finite to small separation (9j HO]. In line with these, it 
should be mentioned that for Neumann boundary conditions the numerical results are 
less reliable and a satisfactory agreement with the analytical results could not be reached 
so far. 

In the present paper we consider analytically a further limiting case, namely a 
cylinder A of fixed radius Ra at finite separation from a second cylinder B (see Fig.l), 
whose radius Rb becomes large, Rb — > oo. In the limiting case we reproduce, of course, 
the result for a cylinder on front of a plane. We derive the general expressions for the 
first two corrections for large Rb and calculate the first order correction explicitly. 

It should be mentioned that the limit of one cylinder becoming large cannot be 
obtained by PFA since the separation between the two cylinders remains finite. As well, 
this limit cannot be obtained by a dipole approximation since arbitrarily high orders of 
the orbital momenta related to the cylinder B contribute. This can be seen below on 
the hand of the approximation of the kernel K^ B (ip, ip'), Eq. (l30|) . 

The paper is organized as follows. In the next section we re-derive the formulas for 
the vacuum energy in the presence of two cylinders. We need these formulas in order 
to introduce the necessary notations adopted to the needs of the third section where we 
consider the limit Rb — > oo. Finally, section 4 contains some discussion and conclusions. 
Throughout this paper we use units with H — c—1. 

2. The basic formulas for two cylinders 

In this section we display the formulas for the vacuum energy of a scalar field obeying 
Dirichlet boundary conditions on two parallel cylinders. The configuration is shown in 
FigfTJ The vacuum energy can be written in the form of a trace of the logarithm, 
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1 f°° 

E=—J duu}Tr]n(6(<p-<p') - M(<p, if/)), (1) 

with 

M(<p,<p') = K^{cp,cp'')K AB {cp''^)K^^')KBA(^,^), (2) 
where integration over doubly occuring angles is assumed and the trace in (pQ) is over 
the angles <p and ip'. In these formulas, the kernels K AB (<p, ip) are the projections of the 
free space propagator 

/dk e ik(x-^) 

which is taken in Fourier representation with respect to the translational invariant 
directions x and x 3 , onto the surfaces of the cylinders. The vectors x and x 1 are 
in the (x\, x 2 )-plane and k = (ki, k 2 ) is the corresponding momentum. For the cylinders 
we use the following parameterizations, 

A : 3*{<p)=( L + RaCOB(P ), B: f» = ( R ^~ l + C ° S ^ V (4) 
Ki ^ J y R A sin <p J ' yy! y R B sin <p I K ' 

so that we have 

K AB ( V ,ij) = A UJ (x A ( V ),x B (ij)) (5) 

and accordingly with A <-» 5. The inverse of a kernel is taken in the sense of an 
operation on the surface of the cylinder, i.e., 

dl f" Kaa(p> f>")K AA {(p" , ip') =6(<p- if') (6) 



must hold. In Eqs. ([T]) and ([2]) the integration over the corresponding angles is assumed. 
The interval for all angular integrations is (p G [0, 2ir]. As well, the trace is to be taken 
in this sense, for instance 

TrM(^')= [*" dtp M(<p,<p). (7) 



In this way, all quantities entering the representation (pQ) of the vacuum energy 
are defined. However, in order to work with explicit formulas one need to change this 
representation by introducing an appropriate basis in which the inverse kernels become 
diagonal. According to the geometry of the considered problem we use the basis 



e ilip 



W) = -7= (8) 

V Z7T 



and the notations 

(l\Kcc\l') = Kcc>u> = J d<pd<p' - % ^^ K cc ,(tp,p'), (9) 

where C and C stand for any of A or 5. In this basis, all quantities become infinite 
dimensional matrices in the indices (I, I') (I = — oo, . . . , oo) and the vacuum energy can 
be rewritten in the form 

i r°° 

E=— dwwTr\R(5 u ,-M ltV ). (10) 
4?r Jo 
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For the needs of the next section it is useful to represent the matrix My/ by 

Mi,v = JC& Ntf (11) 

(the matrix K^ A . IV is diagonal, see below Eq. (fT9l) ). where Nu> can be written in 
coordinate space as 

N Ll , = {l\N{^v')\l') (12) 

with 

N((p, <pT) = J K AB ( V , ^K-Kij, WKbaW, <pT) (13) 

or, in orbital momentum representation, as 

Niji = / ] KaB;1,1"Kbb.j»KbA;1",1' ■ (14) 

i" 

In terms of these quantities, the logarithm in the formula ( ITUl) for the vacuum energy 
can be expanded and one comes to the completely explicit representation 

^ /»00 oo ^ oo oo oo 

E =^ E ••• E M Ml M W2 ...M w . (15) 

s=0 i=— ooZi=— oo Z s =— oo 

As already said, this is a finite expression, i.e., the integration and all summations in 
(IT5"|) do converge for any fixed values of the parameters Ra, Rb and L. 

In the next step we need to remind the known explicit expressions for the matrices 
K~^\.i and Nijr. We start with Kaa-,i,v which is in accordance with and (JS]) given by 

/ e —iltp+il'<p' 
dcpdcp' A^xa^), £4 (<£'))• (16) 
Ltx 

Using the expansion 

e izcoa<p = il M z ) eil(p (I 7 ) 

l=— oo 

of a plane wave into cylindrical ones, Eq. (1161) can be rewritten, 

dk k 



r 

KAA;l,l' = K A A;l Ol,V = 

JO 



Jl(kR; 



2tt 

= I 1 (luR a )K 1 (luR a ). (18) 

The last line follows from Eq. (6.541 1) in [TTJ and is in terms of the modified Bessel 
functions. This kernel is diagonal in the orbital momenta and can be inverted simply 
by 

- K AA, *W = j^rX^RaY (19) 

The corresponding formulas for K^ m v follow by substituting the radius, Ra — » Rb- 

The remaining matrices originate from free space Greens functions having their legs 
on different cylinders, 

Kabw= fdcpdcp' e ^ * A, (£ A (¥>), £ B (y/))- (20) 
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Using the expansion (IT7|) three times and a corresponding formula generalizing 
Eq. (6.541 1) in [11] one comes to 

Kabw = (-l)^M? A )i^KL + R B ))I v {uR B ). (21) 

This formula gives the transition from the cylinder A to the cylinder B if using the 
terminology of the transition formula approach. The reverse formula follows by spatial 
reflection of the plane x\ — (L + Ra)/2, 

K BA ,^ = {-l) l+l ' K AB , V ^. (22) 
In a similar way, with the substitution 

Rb —Rb (23) 
one comes to the formula 

Kabw = {-^I^Ra^i-v^Rb - L))K v {u)R B ). (24) 

which corresponds to the configuration of the cylinder A inside the cylinder B, see Figj2j 
With the above formulas, i.e., by inserting Eq. (1181) and fTSOj) with ([2]) into f|T5|) we 
come to the known formula for the vacuum energy with Dirichlet boundary conditions 
on the two cylinders of radii Ra and R B . As already mentioned, this formula allows for 
a direct numerical evaluation at fixed Ra, Rb and L. 



3. One cylinder becoming large 

In this section we consider the vacuum energy ( ITUl) in the limit Rb — > oo. In the leading 
order we reproduce the corresponding expression for one cylinder in front of a plane. 
The next-to-leading order will then be the main result of this paper. 

We start from representation fTTUT) of the vacuum energy with M\y given by Eq. (1111) 
and the inverse kernel Kaa-,i given by Eq. ffTB]) . All quantities related to the cylinder B are 
contained in N^, Eq. ( 112ft . For the following it is convenient to consider its coordinate 
space representation ({TBI) . The main step for considering the limit Rb — > oo is the 
expansion of the kernel Kg B (tp, ip 1 ). We use its orbital momentum sum representation, 

KbbM') = T(/jR ^ y (25) 

2ir /-^ 1i(uRb)Ki{ujRb) 

l=— oo 

and inserted the analogue to Eq. (jT9l) . Now, if we would take the limit Rb — > oo in the 
last expression, using the asymptotic expansion of the Bessel functions, 

I v {z) ~ e 2 , K v {z) ~ e~\ (26) 

we would get a diverging sum over /. Hence, arbitrarily high I contribute to the 
considered limit. Therefore we have to take the uniform asymptotic expansion of the 
modified Bessel functions |12|, 
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(we do not need the explicit form of the function r](z)). With v — > I and z — > uoRb/1 we 
get 

'<(^>«> - vftW ( l + t2(1 ~^ 1 ~ 5i2> + ' ' ■) (27) 

with 



^p + (uR B y' 

Further, in Eq. (fT3|) we change the variable of integration ip for Z2 according to 

^ = arcsin |l = ^ + i^v + ... ( 28) 



Rb Rb 6 \Rb / 

with 

The range ^2 £ [ — Rb, Rb] will extend to the whole axis in the following. We also make 
the corresponding substitution for the primed quantities in Eq. (lT5j) . In fact, this change 
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of variables is the orthogonal projection of the right half of the circle corresponding 
to the section of the cylinder B onto the axis z\ — 0. It is to be mentioned that the 
left half of that circle does not contribute to the Casimir force in the limit Rb — > oo. 
This statement would not be true for any finite Rb, however we are going to obtain an 
asymptotic expansion. 

With these expansions, the kernel represented by Eq.( |25]) becomes a function of z 2 
and z' 2 i 



OO 

K BB (z, z') = — 2 \A 2 + Ms)* (3°) 

1 fu\Yl - 5^_) + , p ((z2)3 _ (4)3) \ + \ e ^ (z2 _ 4) _ 



2tt 

l=— oo 



x I 1+ - ' 



Rl\ sre 

Finally, in the limit Rb — > oo, the sum over the orbital momenta in (125]) becomes an 
integral after 

1 °° poo j 



t=— oo 



In place of ( 1251) we get 

/oo 
rfg 2T q e^"-^ H(q), (32) 
■oo 



where 
and 



r g = v^T7 (33) 



HUj) 1 + ^| (^^j " + «K*2 - + • • • • (34) 

Now we are going to insert (1521) into iV(</?, </?'), Eq. (lT3"j) , where we make the substitution 
(T28|) . Further we use the momentum space representation 

A^zJ) = r ^ -i- e ^-4)-r fc k-4l ; (35) 



.3 /3x 



00 



2vr 2r fc 



which can be obtained from ([3]) by carrying out the integration over k\ and renaming 
k 2 as fc. For Kab^,^) we use its definition where we substitute x B (ip) — >■ 2 with 
Z2 following from Eq. (1281) and £1 given by 




The result is 

AW) 



Zl = BB |_ 1 + l/l _(^] | = __|_ + .... (36 ) 



dz? f dz' 2 



2 / / / , \ 2 



dk ^H^M-^-Tk^^-z^ r e ifc'(^(v3')-4)- r fe'l a: i l (¥'')-4l 



X 



2vr 2r fc y 2?r 2r^ 

R% J dq 2T q e- q(z2 ~ z '^H(q) (37) 
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(here and in the next formulas all integrations go over the whole axis) . In this expression 
an expansion in 1/Rb can be made. Keeping contributions up to second order we get 

AW) = / 'dzj 'dz'Jf- If- I 'dq H 



with 



27r J 27r J 2T k T k i 

x e iz2(-k+q)+iz' 2 {k'-q)+ikx^( t p)-T k x^( t p)-ik l x^( l p')-T kt x^( } f,') /g^ 



i 1 (T k zl + T k ,z' 2 2 ) | 1 



Rb 2 R 2 B 



(r k zl + r k ,z> 2 2 ) 2 {zl + z' 2 2 ) 



+ .... (39) 



Here we used also xf((p) — Z\ > 0. In (jHSj) the integrations over z 2 and z 2 can be carried 
out delivering delta functions and their derivatives which allow to remove two out of 
the three momentum integrations. 

Let us first consider the zeroth order term, in the expansion which follows 

from fl38l) with H — > 1. It reads 

N^(cp,cp') = [ dq — e ^2^)-^W))-T q {xt(v)+^{v')) _ (40) 

J 

The corresponding quantity in the orbital momentum basis (defined as in Eq.Q) is 

i+i' 

N$ = (-I^'KljRaWuRa) jdq±r (^) ' e~ 2 ^ L , (41) 
where we used 

f * <jV e -ur««4( v )-r q *f( v ) = ^ly^RAe-^L ( ^Vzl) 2 . ( 42 ) 
Jo 2tt \T q + qJ 

The last expression can be derived from (fTTl) . The remaining integration can be carried 
out using Eq. (8.432 1) in [11] and results in 

N $ = (-l) lM 'li(uR A )K l+v (2uL)I v {uR A ), (43) 

in agreement with the corresponding formulas for a cylinder in front of a plane, see An,m', 
Eq.(A7), in [5] or the corresponding formulas in j3j. Being inserted together with ffT9l) 
into ffTUl) with account for ffTTl) the leading order in the limit Rb — > oo reproduces just 
the energy for a cylinder in front of a plane. 

Now we consider the first next-to-leading order. It results from the 1 /-Re- 
contribution in H, Eq. (l39l . and its contribution to N(<p,(p'), Eq. (l38l) . is 

N^Up.ip') = — I dq ( — j^^^A fl e -^M-r^w) . (44) 
2Rb J \dq J \dq J 

In calculating this expression from (138]) we represented z\ by (d/dk)(d/dq) and 
integrated by parts each derivative. The contribution from z' 2 appears to be the same. 
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Now we calculate the corresponding expression in orbital momentum basis. The 
angular integrations can be carried out as before using (|42p . The corresponding integral 
in the angular momentum representation can be written in the form 

i> 



^(-iyi lt ( u R A )[^-JL) e -™<\. (45) 



| ( _ 1)I%MW (^ 

Simplifying this expression we obtain for N^J finally 

N$ = ^ I^RMuRa) 

ZK b 



i+i' 

x ldq±(qL + I) (qL + I') (Ei^J ' e~ 2 ^. (46) 



Here the integration over q cannot be carried out as easy as in Eq. (|4ip and we keep it 
as is. 

In this way we obtained the expansion 

N ltl , = AT$ + AT$ + . . . , (47) 

where N^] , Eq. (j4"3l) . is independent from Rb and N$ , Eq.()45l). is of order 1/Rb- The 
dots denote the contributions of higher orders. Now we insert this expansion into Eq. (jTOl) 
using (JIB and (jj7j), 

E=±- TduuTrln (s u , - K A \^ (jV$ + JVg? + ...)), (48) 
and expand the logarithm, 

£ = — jf du u Tr In (\, - l^X' ) 

- - / u Tr - /qi.^?) K- A \,N£}„ + . . . , 

*/ 

= #(o) + #(i) + . . . . (49) 

The leading order, E^ -* is the energy for a cylinder in front of a plane and E^- 1 ' is the 
correction of order 1/Rb- The latter can be rewritten in the form 

1 POO . \ ~1 

E {1) = -£-\ duuTr (S^Kaai - jV$J Nfi„. (50) 

«/ 

Eq. (!50l) is the final step in the calculation of the l/i^-correction to the vacuum energy 
for the radius of the cylinder B becoming large. It can be calculated numerically since 
all sums and integrations entering do converge. 

In order to represent the result in a more instructive way we represent the energy 
in terms of dimensionless functions, 
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where 

d = L-R A (52) 

is the separation between the two cylinders. For Rb — > oo we rewrite the last line in 
Eq. (1491) in the form 

where E^\d/R A ) is a dimensionless function describing the case of the cylinder A in 
front of a plane and E^'(d/ 'Ra) describes the first correction for large Rb- Further we 
rewrite Eq.( l53l) in the form 

where AE(d/R A ) = E {1 \d/ R A ) / E^\d/ R A ) is the relative correction. 

The behavior of the function E^\d/R A ) is well known. For large argument, i.e., 
for d » R A , it describes a small cylinder (or a cylinder at large separation) in front of 
a plane. It has a logarithmic behavior, 

#o) (*) I?_ (55) 

which is due to the logarithmic behavior of the two-dimensional Greens function ([3]). 
For small argument, i.e., for d << R A , its behavior follows from PFA, 



\RaJ 1920^ V d \ 3QR A J' K J 

where we also included the first correction beyond PFA [!5j. Numerical evaluations of 
this function can be found in [3j [TO] . 

The function AE(d/R A ) can be calculated in a similar way. We start with large 
arguments and consider first the function E^(d/R A ). We expand for small Ra using 

h (ujR a ) K t (coR A ) = 5i, f - In ^ - 7 + . . . ) (57) 



and get from ( fl5 



1 r<x> rP'T^ 

N}$ = koSv fl ^- d q q -^re-^ + .... (58) 



oo 



Here we took into account that we get in this approximation only the contribution from 
I — I' — in A^y, Eq. pBI) . Because of the explicit factor Ra in _D n -(g), Eq. fHBl . only 
i = I' — did contribute. In this way we come to 



With the substitutions 



q = uj sinh 0, = — , (60) 

zL 
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and in leading order for small Ra and with (1521) we get 

Ra) 64irln#i„ J_„ cosh 9 



Ra 

The integrations can be carried out resulting in 

EF> (-*- )~ - — (62) 

UJ 487rln^ ^ 

Together with ( 155]) and ( 1541) we get for the relative correction at large separation 

**(£)-->+ <>$)• 

which is the limiting value for this function for large argument. 

Now we consider the opposite limit of small argument. It corresponds to a large 
cylinder A in close separation from the larger cylinder B. In this case the energy can be 
calculated from PFA. In general, if both cylinders are large, the energy Ea,b is given in 
PFA by 

71-3 1 [r 

Ea - b = -^M^1' (64) 

where 

R = -^*§- (65) 
Ra + Rb 

is the geometric mean of the two radii. It can be expanded for Rb » Ra, 

R = R a -^ + .... (66) 
Kb 

So we get for E, Eq. (}4"9"]) . for d « Ra « Rb, the following approximation, 



F ^ 1 ?± fl 4- - — - l?± \ (RTl 

1920^d 2 V d V S6R A 2R b + ---)' [ ° n 

where we added also the correction beyond PFA which may be of the same order as 
the correction for large Rb- Comparing (|54|) with (IBTl) we infer the behavior for small 
argument, 

a K£h-^t +o(i) - (68) 



As already mentioned, this function can be evaluated numerically using Eq.( 150j) . This 
evaluation was carried out with a truncation of the orbital momenta sums which ensures 
a precision of at least 2 digits. Surprisingly, the numerical results can be described 
remarkably well by the fit 

fi , / d \ Ra ( Ra \ ( Ra ^ 



AE ™ ^—j = 0.359 + 0.514 -f - 0.002 j +0.0001 . (69) 

In Figj3]it is shown how the function dj 'RA^-Eidj Ra) approaches the value 1/2 which 
follows from the PFA, Eq.flEBJ. 
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Figure 3. The function d/R A AE(d/R A ) and its PFA-limit 1/2 



Now we consider the configuration of a cylinder A inside the cylinder B, see Fig. 2. 
The calculation goes in close parallel to the former case and we indicate only the 
necessary changes. First of all, we have to define the coordinates parameterizing the 
cylinder B now. These were given by (0J. Whereas x A does not change, we have for x B 
now 



which appears after a reflection on the plane (x\ = 0). 

The expression fl25|) for K^(ip, ip') does not change. We can keep the substitution 
fl28|) and fl29|) . whereas in place of fl36|) we have now 



The first term of the expansion has the opposite sign as compared to ( 1361) . In fact, this 
is the only change we have to account for in the subsequent formulas. In these we have 
first to consider K^ B [z^ z'), Eq.(l30l). It remains unchanged. Next is N(ip,ip'), Eq.(|371). 
Here the variables Z\ and z[ appear only in the exponential together with and T^/. 
In making in (|37|) the expansion in 1/Rb, the sign change appears in H, Eq. (l39j) . just in 
the first order contribution. All remaining calculations go in the same way as before. In 
this way the changed sign can be traced until the final formula for the energy, Eq. fl63|) . 
which in the case of an inscribed cylinder reads 



In this case the energy is increased, again in correspondence with the expectations. 




(70) 




(71) 




(72) 
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4. Conclusions 

In the forgoing section we considered the vacuum energy of a scalar field obeying 
Dirichlet boundary conditions on two cylinders and calculated the asymptotic expansion 
of this energy for one of the cylinders becoming large. We have shown how to construct 
this expansion and wrote down the first two orders in general form. As a particular 
example we considered the first order in the special case when the separations between 
the cylinders becomes large, Eq. (l63l . The other case, when the separation becomes 
small is covered by PFA and resulted in Eq. fl68l) . 

The asymptotic expansion for large radius Rb involves arbitrarily high orbital 
momenta for the cylinder B. This is similar to the expansion for small separation, 
but in detail, of course, different. It should be mentioned that the limit of one cylinder 
becoming large cannot be obtained from PFA unlike the case of small separation. In this 
sense it is an independent calculation. However, it should be related to a perturbative 
expansion which emerges if considering the large cylinder as small deviation from a 
plane. For consistency reasons, it would be interesting to check this. 

V.N. was supported by the Swedish Research Council (Vetenskapsradet), grant 621- 
2006-3046. 
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